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Abstract. Let / be a degree D univariate polynomial with real coefRcients and exactly m 
■ monomial terms. We show that in the special case to = 3 we can approximate within e all the 

\ roots of / in the interval [0, R] using just 0{\og{D) log [D log -p)) arithmetic operations. In 

particular, we can count the number of roots in any bounded interval using just 0(log^ D) 
arithmetic operations. Our speed-ups are significant and near-optimal: The asymptotically 
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' sharpest previous complexity upper bounds for both problems were super-linear in D, while 

our algorithm has complexity close to the respective complexity lower bounds. We also dis- 
cuss conditions under which our algorithms can be extended to general m, and a connection 
to a real analogue of Smale's 17— Problem. 

< 

ri . 1. Introduction 

^ ■ 

^ ' Real-solving — the study of solving systems of polynomial equations over the real numbers 

— occupies a curious position within computational algebraic geometry. From the point of 
view of computational complexity, classical algebraic geometry has left real-solving almost 
^ ! completely untouched. For example, rigorous lower bounds for the arithmetic complexity of 
\ finding approximations to the roots of polynomial systems didn't appear until the work of 
■ Renegar in the late 1980's |Ren87| IRen89j . and finding optimal bounds continues to be the 
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subject of much active research |BMS1'97[ IB(LSS981 IMP99L |Roj99bl KiLSOlL IMPH,n3j . 



O '. As for counting real roots, there is a beautiful result that one can bound their number 

independently of the degrees of the underlying polynomials: Askold Khovanski proved an 
explicit upper bound singly exponential in the number of variables and the total number 
of monomial terms |Kho91j . While this fewnomial bound is far from optimal in higher 
dimensions |LRW03j . it is significantly smaller than the number of complex roots when the 
underlying polynomials have sufficiently high degree. (Sparse polynomials are sometimes 
also known as lacunary polynomials and, over M, are a special case of fewnomials — a more 
^ . general class of analytic functions of parametrized complexity |Kho91j .) 

We are then naturally lead to suspect that a similar improvement is possible for the harder 
problem of approximating the real roots. So can one solve sparse polynomial systems over 
the real numbers significantly faster than via the usual algorithms based on complex algebraic 
geometry? The existence of general speed-ups of this nature is still an open problem, even 
in the univariate case: For example, until the present paper, it was still unknown whether 
the real roots of a univariate trinomial of degree D could be approximated within a number 
of arithmetic operations sub-linear in D [MP99i |Roj99bt K^LSOll IMFRHSj . 



We now answer this last question affirmatively as follows: Let / C C be any subset, 
/ G C[xi] a degree D polynomial, and suppose Ci, . . . , Cm are all the distinct roots of / in the 
region /. By e-approximating the roots of / in / we will mean finding complex numbers 
zi, . . . ^ Zm&I such that for any i G {1, . . . , m}, we have \Q — Zj\<e for some j G {1, . . . , m}. 
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(So ^-approximating in I implies a correct root count in I as well.) When ICR and / G 
we will further stipulate that the Zi all lie in M. 

Theorem 1. Let R,e>0 and suppose /gM[xi] \ {0} has degree D and at most 3 monomial 
terms. Then we can e-approximate all the roots of f in the closed interval [0, R] using 
just 0(log(D) log log j)) arithmetic operations. In particular, we can count exactly the 
number of roots of f in any bounded interval using just 0{log^ D) arithmetic operations. 

As is standard, we count arithmetic operations as field operations in the field over Q gen- 
erated by the coefficients of /. Note also that our underlying algorithm handles degenerate 
roots (i.e., roots of multiplicity >1) with no difficulty. 

Remark 1. Throughout this paper, all 0-constants and Q-constants are absolute and effec- 
tively computable, o 

Remark 2. Our speed-ups are significant: The current asymptotically sharpest (sequential 
worst-case) arithmetic complexity upper bound for e- approximating the roots of a general de- 
gree D univariate polynomial in the open disc {z \ \z \<R} IS 0[D\og^{D)\og\og^) BPSl 
INR96j . The analogous upper bound for counting real roots in an interval is 0(p\o^{p) log log/)), 
via the technique of Sylvester-Habicht sequences |Roy96t ILMOl' j. In particular, no 
sharper upper bounds were known before for the problems considered in Theorem^ above, o 

Remarks. Our speed-ups are also near-optimal: For any fixed D, R> 2, our e- approximation 
algorithm matches ( up to an asymptotically constant multiple ) the VL (log log -) arithmetic 
complexity lower bound known for e- approximating the roots of x"^ — N where 1 < N <2 
|BMST97] . As for counting the roots, there do not appear to be any explicit lower bounds 
known for the arithmetic complexity of counting the real roots of m-nomials. However, one 
should note that the arithmetic complexity of just evaluating a degree D monomial is Q{logD) 
in the worst case |dMS96| IMor97j . o 

Our algorithms are based on an earlier hybrid algorithm of the second author which 
combines bisection and Newton iteration, a new observation on the Sturm sequences of 
trinomials (Theorem^lof Section IT!^ . and some analytic estimates (Theorem^lof Sectional). 
We in fact give a more general algorithm, in Theorem IHl of Section |2l that applies to certain 
univariate m-nomials. However, for m>4, there are two main obstructions to showing that 
our more general algorithm has complexity sub-linear in D. We refer to these obstructions 
as Problems A and B, and describe the first now. 

Problem A. Is there an absolute contant k such that one can count exactly the number 
of real roots (in any input interval) of arbitrary m-nomials of degree D using, say, just 
0{\og'^"^ D) arithmetic operations? o 

One may note that we have hedged our bets in Problem A by letting the complexity bound 
increase exponentially in m. This is motivated by the NP-hardness of the multivariate 
analogue of Problem A. 

Proposition 1. Suppose f G Z[xi, . . . is an m-nomial and we measure the size of f 
as the total number of bits necessary to write the binary expansions of all the coefficients 
and exponents of f . Then it is NP-hard (in the classical Turing sense, using the preceding 
notion of size) to decide whether f has a real root or not. ■ 
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(The proof follows immediately from a standard reduction of the special case of n-variate 
quartic polynomials to SAT |BSS89j .) We present some examples in Section [1.21 revealing 
that the classical approach of Sturm sequences will most likely need to be abandonded for 
m > 3. Nevertheless, speed-ups for m fixed still appear possible and are of considerable 
practical interest for large D. 

Noting that ^-approximation extends naturally to polynomial systems (by ^-approximating 
each coordinate separately), a consequence of our univariate trinomial algorithm is the fol- 
lowing result which may be of use for solving general pairs of bivariate trinomials. 

Theorem 2. Suppose f G ]R[a;i,X2] has exactly 3 monomial terms, D is the degree of f, 
and5>0 is any constant. Then, using just O (log{D) log (^D log j)) arithmetic operations 
and O (log^''''^D) bit operations, we can e-approximate all isolated inflection points, vertical 
tangents, and singular points of the curve {(a;i, ^2) G [0, i?]^ | /(xi, 0:2) =0}. 

The asymptotically sharpest (sequential worst-case) arithmetic complexity upper bound for 
^-approximating the roots of a system of n polynomial equations in n variables are at best 
polynomial in log log — and a quantity sometimes attaining (see, e.g., |MP99| Roj99b 



IGLSOU IMPR03j ). The importance of Theorem |21 lies in that it is a first step toward a 
practical two-dimensional analogue of bisection. The latter in turn is a first step toward 
generalizing our algorithms here to higher dimensions. 

We leave the bit complexity of our algorithms for a future paper. In particular, unless 
stated otherwise, our underlying model of computation here is the BSS model (with inequal- 
ity) over M [BCSS98j . Those unfamiliar with this model can simply think of such a machine 
as an ordinary Turing machine augmented with registers that allow arithmetic and inequal- 
ities with real numbers as well as bits. However, let us touch upon a deeper question: The 
connection of our work to Smale's notion of approximate roots (see, e.g., |Sma86t lBCSS98j 
and Section 121 below). 

Definition 1. Suppose fi, ■ ■ ■ , fn^ M[xi, . . . ,Xn], F := (/i, . . . , /„), and the total number of 
distinct exponent vectors in the monomial term expansions 0//1, . . . , /„ is m. We then call 
F a real m-sparse n x n polynomial system, o 

Definition 2. Let B C<C"' be any open ball and F : B — > C" an analytic function. Also, 
given any zqEB, let us define the sequence (-2j)^o ^i-=Zi-i ~ J^c(F)|~^^F(zj_i) for all 
i&N. Then, if there is a root (^B of F such that \zi+i — (\ <8(|)^'|2;o — CI f(^f^ (ill i (i.e., if 
Newton iteration for F , starting at zq, converges quadratically to a root of F), we call Zq an 
approximate root of F with associated root (. o 

Real Analogue of Smale's 17— Problem. For fixed n, can all real roots of an m-sparse 
n X n polynomial system be found approximately, on the average, in polynomial time with a 
uniform algorithm? More precisely, let F be an m-sparse nxn polynomial system with max- 
imal exponent D and coefficients that are, say, independent standard real Gaussian random 
variables. Is there a uniform algorithm that finds a set of approximate roots close to all the 
real roots of F, with average-case arithmetic complexity 0{{mlogDy) for some constant v 
depending only on n? o 

The polynomial system ((xi — l)(xi — 2), ... , (x„ — l)(x„ — 2)) clearly shows that fixing n 
is necessary in our real analogue above. The appellation uniform merely emphasizes that 
there be a single algorithm, with all steps explicit and constructive, which works for all 
inputs |B(XSS98j . 
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The original statement of Smale's 17— Problem |Sma98j (see also |SmaOO| Pg. 287]) differs 
from our analogue above as follows: (1) n is allowed to vary (so one seeks an absolute 
constant z/), (2) one averages over choices of complex coefficients, and (3) one instead asks for 
a single complex approximate root. Smale also left the underlying probability distribution 
unspecified. We observe that a positive answer to the following variant of Problem A would 
be quite useful in the direction of our real analogue: 

Stochastic Version of Problem A: Suppose f is a univariate m-nomial of degree < D with 
coefficients that are real standard Gaussian random variables. Is there a uniform algorithm 
that counts exactly the number of roots of f (in any input interval), with average-case 
arithmetic complexity 0{{mlogD)'^) for some absolute constant k? o 

Smale's original 17— Problem remains unsolved, although a partial affirmative answer (an 
algorithm containing a non-constructive step which may be called many times) was found 
by Shub and Smale in the mid-1990's ^Ml llj(HS98l . 

1.1. Earlier Work on Approximating Roots. 

When can we solve a polynomial system in time polylogarithmic in the degree of the under- 
lying complex algebraic set? An affirmative answer is trivial for the special case of a single 
polynomial with < 1 monomial term, and was known at least since the mid-1970's for the 
(univariate) binomial case, e.g., |Bre76j and |Ye94| Sec. 4]. On the other hand, little seems to 
be known about the case of 3 or more monomial terms: To the best of the authors' knowledge, 
the only result close in spirit to Theorem^ is a result of Daniel Richardson |Ric93j implying 



that the arithmetic complexity of counting the number of real roots of ci + C2x'^ + c^x^ (with 
0<d<D) is polynomial in d and logD. 

Not much else seems to be known about the intrinsic complexity of real solving, or even 
real root counting, for univariate polynomials with 4 or more monomial terms. One known 
speed-up over the usual univariate algorithms over C is |RojOOa[ Main Theorem 1.2] which 
gives an arithmetic complexity bound of 0(mlog(-D) log (D log j)) for e-approximating the 
roots of 

CiX"' +■■■ + CkX"" - Cfc+iX"'=+^ CmX""^, 

when < Oi < ■ ■ ■ < = -D are integers and all the Cj are positive real numbers. 

It is also not difficult via the results above to construct various systems of multivariate 
trinomials which admit super-fast solving in the sense of Theorem ^ as well. For example, if 
one has /i eM[a:i], /2 gM[xi, X2], . . . /„ gM[xi, . . . , x„], and all the are trinomials of degree 
< D, then to solve F := (/i, . . . , /„) one can simply solve /i first and then recursively solve 
the resulting smaller system. Letting A/f denote the number of roots in the nonnegative 
orthant of such an F, it is easily checked that A/> < 2" (see, e.g., jLHWn3[ Thm. 3, Part 
(c)]). We can then easily derive, via Theorem Q an arithmetic complexity upper bound of 
O (A/> log^ (D) log log f ) for the more general problem where we instead ^-approximate all 
the roots of F in the orthant-wedge 

Vr^:={(xi,...,x„)GM" \xl + --- + xl<R^ , Xi>0 for all i] . 

Extending these results to general trinomial systems, not to mention general sparse systems, 
remains an open problem. Nevertheless, for a general binomial system F with exactly T><oo 
roots in C" (counting multiplicities), one can e-approximate all its roots in using just 
©((logP) (B^log^(n) + loglog j)) arithmetic operations and Oln^B"^^^ \og^^^ n) bit operations, 
where B is the total number of bits needed to write down the exponents of F and 5 > is 
an arbitrary constant |RojOOa[ Main Theorem 1.3]. (The bound in |RojOOa| is written in a 
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slightly different manner and we have here talcen the liberty of improving the bit complexity 
portion by employing a recent result of van der Kallen on computing the Hermite normal 
form of an integral matrix |vdKnnj .) 

1.2. A New Algebraic Observation and Earlier Work on Counting Roots. 

As for merely counting the real roots, there are still surprising gaps in our knowledge. 
For instance, the main general algebraic techniques for real root counting — Sturm-Habicht 
sequences |Htu351 IHab481 |Roy96| ILMnij and Hermite's Method |Her56[ IPH,S93[ |Roy96| — 



originated in the 19— century and are still hard to improve upon. 

Definition 3. For any sequence of real numbers s := {si, . . . , Sk) , the number of sign 
alternations of s, Ng, is the number of pairs with l<j < j'<k, SjSjr<0, and Sj = 

for all iE{j + I, . . . ,j' — 1}. Also, for any polynomial /gM[xi], define the sequence (pt)^o 
where po := f , Pi-=f', Pi+2-=qiPi+i—Pi, Qi is the quotient of and px is the last element 
of the sequence not identically equal to 0. We then call {po, . . . ,pk) the Sturm sequence 
off. o 

Lemma 1. |Stu35l |Roy96| Following the notation of Definition]^ let a, 6 G M with a<b, 
A:= {po{a), . . . ,Pk{(i)), and B := {po{b), . . . ,pK{b)). Then the number of distinct roots of f in 
the open interval (a, b) is exactly — Nb ■ In particular, /ai < ■ ■ ■ < and Ci, . . . , G M/ 
=^ Cix""^ + ■ ■ ■ + CmX"""" has no more than N(^ci,...,cm) (<rn — 1) positive roots. ■ 

Remark 4. The latter part of the above lemma is Descartes' Rule of Signs which dates 
back to the famous philosopher's 1637 book La Geometric (see also |,SL54. Pg. 160]j. o 

The following observation on the Sturm sequence of a trinomial, which we use in proving 
Theorem^ may be of independent interest. 

Theorem 3. Following the notation above, suppose f has at most 3 monomial terms. Then 
i^r< 3[log2 -D] +2 and {pi, . . . ,pk) consists solely of binomials, monomials, and/ or constants. 
In particular, the entire Sturm sequence of f can be evaluated at any real number using just 
0{\og^ D) arithmetic operations. 

Extending our last theorem to polynomials with more monomials appears unlikely: First, 
note that the quotient of the division of two binomials can be quite non-sparse, e.g., ^ = 
+ ■ ■ ■ + a; + 1. So the expansion of any intermediate quotients must be avoided. Further- 
more, the tetranomial case already gives some indication that Sturm sequences may in fact 
have to be completely abandoned: the following example shows that, for all D>2, the fourth 
element of the Sturm sequence of a degree 2D tetranomial can already he a {D + l)-nomial. 

Example 1. Consider the tetranomial po{x) := + x^~^^ + a;'^ + 1 and its derivative 
Pi '■=p'q = 2Dx'^^~^ + + l)x^ + Dx^~^ . The resulting Sturm sequence then continues with 
P2 = """^^3;^"^^ — — 1 and from here it is easy to see (by a writing a simple recursion 
for the resulting long division) that the quotient of ^ is a polynomial of degree D — 2 
with exactly D — 1 monomial terms. Thus, ps :=q3P2 — Pi has degree D and at least D + 1 
monomial terms, o 

We are also willing to conjecture that the maximal length of the Sturm sequence of a degree 
2D tetranomial is Q[D). (Maple experiments have verified this up to D = 150.) Nevertheless, 
while the behavior of tetranomials is thus more complicated, this example need not rule out 
a more clever method to circumvent these difficulties. 
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1.3. A Connection to Discriminants. 

If one insists on relying on Sturm sequences then one is naturally lead to the univariate 
sparse discriminant. Briefly, given an ra-variate m-nomial / with indeterminate coeflicient 
vector C and set of exponents A, its sparse discriminant (or A-discriminant) , A^(/), is 
the unique (up to sign) irreducible polynomial in Z[C] \ {0} of lowest degree which vanishes 
whenever C is specialized so that / has a root in (C\{0})" in common with . . . , 
|GKZ94j . It is not hard to see that if (a) n = l, (b) the coefficients of / are all constants, 
and (c) A is exactly the set of exponent vectors of /, then we have A^(/) = <^==^ = 1 in 
the Sturm sequence {po, . . . ,pk) of /. So the following is clear. 

Proposition 2. Suppose /gM[xi] has set of exponent vectors A and one can compute the 
signs of the Sturm sequence of f evaluated at anyr&'R using just T{A) arithmetic operations. 
Then one can decide the vanishing of AA^f) using just T{A) arithmetic operations. ■ 

It would be quite enlightening to understand the converse of Proposition |21 and its possible 
obstructions. 

Perhaps more than coincidentally — recalling Theorem Q — ^a(/) can be computed in 
polynomial time when / is a trinomial: up to sign, the formula is simply 

when 02 and 03 are relative prime and < 02 < a3 (see jC;KZ941 Pg. 406] and Remark [7| 
of Section El below). Whether polynomial time complexity persists for general m-nomials 
is already non-trivial for the tetranomial case since the underlying discriminants become 
much more complicated. For example, up to sign, Ajo,6,io,3i}(ci + C2X^ + c^x^^ + c^^x^^) is the 
following homogeneous 21-nomial of degree 35: 

_2i032i5i072i^4^3i _ 21332455721^2^5^28 _ 214327721^10^25 ^ 28315567" • 19 ■ 3l6ll731clcicf c2 
-2*3!^ ■ 5 ■ 7^^ ■ 19 ■ 31"* ■ 4931 ■ 11924843839cf c|c§Oc| - 2^'i^^^^'7^''^?,l^^c\^C2cl^c\ - 2*3^6512713 . ^^7 . ^^2 . 3629537cf c^^c^'^cl 
-263^ ■ 5 • 7" ■ 29 ■ 31" ■ 6361 ■ UniG'icf' clef c\ - 2^^i^'^5'^'^7'^'-^ cic\^clAcl - 2^?,^b^'^7'^?,l^ ■ 857 ■ 295269701cf c^^cj^cl 
-2'?!^ ■ 5 ■ 7^31^8 . 3327253cf ct4°c| - 263*5227'i3l''l60730667473cf c|ct + 2H^-o^'^l^ ■ 11 ■ 29 ■ 31^^ . 233 . 15590^0^44 
-265^^31'^ ■ 248552699041cf circlet - 2^3*'57'' • 11 ■ 29 ■ 31^''cf°c|c|c| - 2^ • 35^^31" ■ 113 ■ 317 • 461c|[^4*c|c| 
-2^335423^3 . 41 . 367c2c26cic| - 26 ■ 3 • . 19 ■ 31^3^8 4cic| + 24345^23112^9^19^^^6 _ 3181^25^10 _ 2636550ciict 

whose largest coefficient has 47 digits. So we pose the following problem: 

Univariate Discriminant Complexity Problem. Given any finite subset A C Z, define 
its size to be the number of bits needed to write the binary expansions of all the points of A. 
Can one decide the vanishing of A-discriminants (for, say, specializations of the coefficients 
in C) using a number of arithmetic operations polynomial in the size of A? o 

This highlights an embarassing gap in what is known about the complexity of discriminant 
computation: it is known that sparse discriminants for bivariate polynomials can be com- 
puted in polynomial time on a Turing machine (resp. BSS machine over C) only if NP C P 
(resp. NP C BPP) [RojOOat Main Theorem 1.4], but no such hardness result is known for 
the univariate case. The preceding inclusions (see |Pap95| for a beautiful introduction to 
complexity theory) are currently considered quite unlikely. On the positive side, it is known 
that y4-discriminants (for arbitrary finite A C Z") can be decided in polynomial time when 
A has n + 1 or fewer points |GKZ94t Prop. 1.8, pp. 274-275]. Also, the vanishing of A- 
discriminants can be decided within the complexity class under a number-theoretic 

hypothesis strictly weaker than the Generalized Riemann Hypothesis |Roj04a| . 
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2. Speed-Ups Through a Variant of o-theory 

Checking whether a given point is an approximate root of a given polynomial can be done 
quite efficiently, thanks to the seminal work of Smale in |Sma86j . Let us now formalize a 
refined version of this fact. 



Definition 4. For any analytic function f : C — > C, let 7(/, x): = sup 

k>2 



k\f'{x) 



1 



Remark 5. It is worth noting that l/7(/, xq) is a lower hound for the radius of convergence 
of the Taylor series of f about xq, so 7(/, xq) is finite whenever f is nonsingular at xq 
|B(;SS98[ Prop. 6, Pg. 167]. o 

Recall that a function g : U — > M, defined on some connected domain f/ C M is convex iff 
g{\x + (1 - \)y)<\g{x) + (1 - X)g{y) for all x,yeU and Ae[0, 1]. o 

Theorem 4. |Ye941 Thm. 2] Suppose f is convex and analytic on the open interval {0,R). 
Also let rE{0,R) be any root of f and suppose that there is an a>0 such that X'j{f,x)<a, 

for all X in the closed interval , . Then for all zE (0, R): 

(1) /// is monotonically decreasing in [z, r] and r G \z, (1 + ^)z\ then z is an approximate 
root of f . 

(2) /// is monotonically increasing in [r, z] andrE [(1 — g^)^;, -z] then z is an approximate 
root of f . ■ 

Proposition 3. |Ye94| Prop. 1] Suppose fi, f2 : U — > M satisfy the hypotheses of Theorem 
^ with ai and respectively playing the role of a, and that fx and f2 are either both 
increasing or both decreasing. Then fi + f2 also satisfies the hypotheses of Theorem^ with 
a = ma.x{ai,a2}- B 

Proposition 4. For any D G R\{0} and c G M, the function f : [0, oo) — > R defined by 
X ^ x^ — c satisfies X7(/, x) = | | for all x>0. ■ 

The second author observed the special case D gZ\{0} earlier in |Ye941 Example 1]. The 
proof their applies to aribtrary real D as well. 

An additional technical result we will need is an analytic estimate which globalizes Theo- 
rem 0] But ffist let us observe an oscillation property we will need for our main algorithm. 

Definition 5. For any D,mGN define T{D,m) to be the family of t-nomials f of degree 
<D with t<m. For any f &J-'{D, m)\{0}, let Mo(/) < ■ ■ ■ <'U7v(/)(/) be the ordered sequence 
comprised ofO, the roots of f in (0, oo), and oo. We then say that f is dampened iff for 
all iE{0, . . . , N{f) — 1}, /" has at most one root in {ui{f),Ui+i{f)) . o 

We of course have J-'{Di, mi) ^J^{D2, 7712) whenever Di < D2 and mi <m2. Rolle's Theorem 
applied to /' tells us that («»(/), Mj+i(/)) (for i G {1, . . . , A^(/) — 2}) always contains at least 
one root of /". It is then easy to show via Descartes' Rule and a routine calculation that 
m-nomials are always dampened when m<4. 

The invariant defined below generalizes the local quantity maximized in Theorem H] above, 
and thus helps enforce the accelerated convergence of Newton's method in certain cases of 
interest. 

Definition 6. For any D,m>2, any dampened f G J^{D, m), and i G {1, . . . , N{f) — 2}, 
define Vi{f) to be the unique root of f" in {ui{f),Ui+i{f)) , f-i(/) :=0, and VN{f){f) :=oo. 
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Also define fo(/) to be Ui{f) or the unique root of f" in {uo{f),Ui{f)) , according as f" is 
non-vanishing in («o(/), ■"!(/)) or not. Also define fAr(/)-i(/) to be oo or the unique root 
of f" in {uN{f)-i{f),UN{f){f)), according as f" is non-vanishing in (uiv(/)-i(/), Uiv(/)(/)) or 
not. Finally, define a{D,m) to be 

sup sup sup {\x - Ui{f)\j{f,x)} . o 

f&T{D,m) ie{0,...,Af(/)-l} a;e(t;,_i(/),i>,(/))\K(/)} 

/ dampened contains a root of f 

/(0)^0 and /'(0)=0 

While a{D,m) may be hard to compute exactly, we can at least bound it above and below 
explicitly in some cases of interest. 

Theorem 5. For any D,m>2 we have a{D,m) <a{D,m + 1) . In particular, 

^^^ = a{D, 2) < a{D, 3)<{D- 1){D - 2)/2. 

One application of our a-invariant is the following generalization of Theorem Q 

Definition 7. For any m-nomial f , let 6{f) be the smallest exponent appearing in f . Let us 
then define the operators S : J^{D, m) — > ^{D, m), Li : T{D, m) — > J^{D — 1, m — 1), and 
L, ; ^(D, m) ^(D -2, m- 1) by : = x-^(/)/, L,if) := £S{f), and L^if) := £Sif) . 

Theorem 6. Following the notation of Theorem^and Definition\^ suppose f &J^{D, m) and 
that all the polynomials {{So Lg^ o ■ ■ ■ o Le^){f) | A; G {0, . . . , m — 1}, Cj G {1, 2}} are damp- 
ened. Then one can find e- approximations to all the roots of f in [0,R] using just 

O (^2"^ jlogp) log (^a{D, m) log + m^lC{D, m) 

arithmetic operations, where }C{D, m) is any upper bound on the worst-case arithmetic com- 
plexity of counting the roots of an arbitrary f &J^{D,m) in an arbitrary interval. In partic- 
ular, if f ^J^{D,4:) then the polynomials 

{5(/),^(L,,(/)),^(L,,(L,,(/))),^(L,,(L,,(Le3(/)))) I e,G{l,2}} 
are all dampened, for any DeN. 

It is easily checked (recalling our observations after Definitional) that any m-nomial with 
m<4 satisfies the hypothesis of our complexity bound above. Theorem IHl thus opens up the 
possibility of super-fast m-nomial solving when m > 3 is fixed. In particular, good upper 
bounds on }C{D,m) and a{D,m) appear to be a fundamental first step. 

Problem B. Is there an absolute contant k' such that a{D,m)=o(\og'^™' ? o 

An important related problem, in the spirit of |MR04j . is how to bound from below the 
probability that a random m-nomial of degree D (with m > 4 fixed) is dampened. 

We state our underlying algorithms in the next section. We then prove Theorems ^ IHl 
and El in Section HJ We conclude with the proofs of Theorems |21 and El in Sections El and IHl 
respectively. 

3. The Algorithm and Subroutines 

The central algorithm we use to prove Theorems Q and IHl MNOMIALSOLVE, is detailed below. 
Succinctly, the key idea behind MNOMIALSOLVE is to subdivide the input interval into sub- 
intervals on which ±/ is convex and monotonic, along with some additional sub-intervals on 
which / is less well-behaved. 
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Figure 1 [—R, R] can be expressed as a union of 6 intervals on which it/ is convex and monotonic, 
and 5 additional intervals of width <2e containing the roots of /'/". o 

Then, using special properties of sparse polynomials, a suitable combination of bisection and 
Newton iteration (detailed in the subroutine HYBRID) yields ^-approximations to all the roots 
in (0, R). (Checking whether and/or R are roots can be done simply by evaluating there.) 
Each approximation is guaranteed to correspond to its own unique root, and the roots in 
intervals on which / behaves badly are certified via a subroutine called FASTERCOUNT. The 
aforementioned subroutines are described shortly after our main algorithm. 



ALGORITHM MNOMIALSOLVE 



INPUT Real numbers R and e with < e < i?, a univariate polynomial f{x) ■.=cix"'^ + • ■ ■ + CmX°'"^ £ 
M[x] with exactly m monomial terms (where < ai < • • • < = D) and (following the nota- 
tion of Definition[7I) all the polynomials {{S o L^.^ o • • • o L^^){f) | /c G {0, . . . , m — 1}, Cj € {1, 2}} 
dampened, and an upper bound a* on a{D,m). 
OUTPUT A (possibly empty) multiset ZCR, such that every root of / in [0, R] is e-approximated by 
a unique Zi € Z, and #Z is exactly the number of roots of / in [0, R\. 

DESCRIPTION 

Step Set Z:=0. 

Step 1 If m = then set Z:=R and STOP. 

Step 2 If m = l and ai>0 then set Z:={0} and STOP. 

Step 3 If m = l and ai = then STOP. 

Step 4 If ai >0 then set Z :={0}, f := f / x""^ , Oj : = aj — ai for all i>2, and then ai :=0. 
Step 5 If am = ^ then set Z:={— 01/02} or 0, according as ciC2<0 or not, then STOP. 
Step 6 If am = '2 then make the change of variables x:=x — 
Step 7 If /(i?) = then append Rio Z. 
Step 8 If m = 2 and f{0)f{R)>0, then STOP. 

Step 9 If m = 2 and f{0)f{R) <0, find an e- approximation z to the unique root of / in (0, R) via 
subroutine HYBRID (with input (e, R,zizf,a*), where the sign is chosen so that / is convex 
on (0,i?)), set Z:=ZU{z}, and STOP. 
Step 10 Otherwise, if 02 > 1, find via algorithm MNOMIALSOLVE an ordered sequence uq< ■ ■ ■ <Uk^ 
comprised of 0, a set of e- approximations (of the correct cardinality) in (0, R) for the roots 
of /' in (0,i2), and R. Then GOTO Step 12. 

Step 11 If 02 = 1 then redefine f{x) :=x"^/(l/x) and £:=min |l, ^, and make the change of 

variables x:=l/x. 
Step 12 Define v^-^ ■.=R. 

Step 13 Using algorithm FASTERCOUNT, decide if /" has a root in [uq + e,ui - e]. If so, let L := 
ui — uq — 2e. Then, evaluating /" at uq -|- e and ui — e, and using algorithm HYBRID with 
input (e, L, ±f, a*) (the sign chosen so that ±/ is convex on {uq + e,ui — e)), define vq to 
be the e-approximation found for the unique root of /" in [uq + e,ui — e]. Otherwise, define 
Vq to be Ui. 



10 



J. MAURICE ROJAS"^ AND YINYU YE+ 



Step 14 Using algorithm FASTERCOUNT, decide if /" has a root in [uk^-i + s,Uk^ — e]. If so, let 
L:=Uk^ - Uk^-i - 2e, g{x) :=/(x + Ufci-i). Then, evaluating /" at Uk^^i + e and - e, 
and using algorithm HYBRID with input (e, L, ±(7, a*) (the sign chosen so that -og is convex 
on + e, Ufcj — s)), define Vk^-i to be the e-approximation found for the unique root 

of /" in [ufci-i + e,tifci — e]- Otherwise, define Vk^-i to be R. 

Step 15 For all i € {1, . . . , /ci — 1}, count via algorithm FASTERCOUNT the number of roots of / 
in {ui — e,Ui + e) n (0, R) closest to Ui and no other root of /'/", and append rrii copies of 
Ui to Z. 

Step 16 For alH E {0, . . . , A;i — 1}, count via algorithm FASTERCOUNT the number of roots of / in 
{vi — e,Vi + e) n (0, R) closest to Vi and no other root of /'/") and append copies of Vi 
to Z. 

Step 17 FOR « G {0, . . . , /fci - 1} DO 
Step 17(a) Evaluate / at + e and — e. If /(uj + e) = (resp. /(uj — e) = 0) then append + e 

(resp. Vi -£)toZ and GOTO Step 17(c). 
Step 17(b) li f{ui + e)f{vi — e)<0 then let Li :=Vi — Ui — 2e and gi{x):= f{x + Ui + e). Then, using 
subroutine HYBRID with input (e, Li, igi, a*) (the sign chosen so that gi is convex on 
(0, Vi — Ui — 2e)), find an e-approximation in [ui + £,Vi — e) to the unique root of / in 
(ui + £,Vi — e), and append this approximation to Z. 
Step 17(c) Evaluate f at Vi + e and Uj+i — e. If f{vi + e) = (resp. /(uj+i — e) = 0) then append 

Vi + e (resp. itj+i — e) to Z and GOTO Step 18. 
Step 17(d) If ^(vi + e)f{ui+i — e) <0 then let Li := Ui+i — Vi — 2e and gi{x) := /(uj+i — £ — x). 

Then, using subroutine HYBRID with input (e, Lj, ztf^j, a*) (the sign chosen so that gi 
is convex on (0, Uj+i — Vi — 2e)), find an e-approximation in (uj -|- e, Uj+i — e) to the 
unique root of / in {vi + e, Uj+i — e), and append this approximation to Z. 
Step 18 END FOR 
Step 19 STOP 

Remark 6. In reality, one should replace the use of subroutine FASTERCOUNT above by the 
best algorithm attaining the complexity bound of}C{D,m) alluded to in Theorem\^ whenever 
m>4. We have omitted this detail above simply to keep the algorithm definite, o 



SUBROUTINE HYBRID (Compare lYeOil Pg. 277]) 



INPUT e,R£R (with 0<e<R), a monotonic analytic function (j) : (0, i?) — > R (with (j}{e)(j){R) <0 and (j) convex 
on (0, R)), and a positive upper bound a* on 27(0, z) valid for all z £ (0, R) . 
OUTPUT An e-approximation of the unique root of (j) in (0,-R), using no more than 0(log (alog evaluations of 
and 0', and 0{\og (a log additional arithmetic operations. 
DESCRIPTION 

Step Define co to be 1 + or — ^— , according as is decreasing or increasing on (0, -R) . Then define Ci :— c^-\ 

for all i > 0, and compute Ci , . . . , cm where M is the first positive integer with cm = Cq > Finally, set 
x:=e and k:=M. 

Step 1 If [(0 is decreasing on (0,-R) and (l){cf,_-^x) > 0) or (<j!> is increasing on (0,-R) and <l>{cf,_^x) < 0)] then set 

a;:=Cj._-^i and k~k — 1 and GOTO Step 1. 
Step 2 Otherwise, if fc>0 then set fc — fc - 1 and GOTO Step 1. 

Step 3 Perform logj (3 + logj ■§) iterations of Newton's method (with x as the starting point), then OUTPUT the 
very last iterate. 

Remark 7. A well-known trick we will use frequently and implicitly is the computation of 
x'^, where x G M and d eZ, within 2\\og2{\d\ + 1)] multiplications using < \\og2{\d\ + 1)] 
intermediate real numbers. Briefly, assuming the binary expansion of d is {ak ■ ■ -00)2 (so oq 
is the I's bit), we simply compute x'^\ via recursive squaring, for all i with ai nonzero, and 
then multiply the resulting numbers together. For instance, 363 = 2^ -f 2^ + 2^ + 2^ + 2 + 1 = 
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(101101011)2 and thus X363 = (. . . (x^) . . . )2 . (. . . (x^) . . . )2 . (. . . (a;2) . . . )2 . ((a;2)2)2 . ^2 . 

8 6 5 

which only requires 17 multiplications and enough memory for 9 intermediate real numbers 
— much faster than the naive 362 multiplications, o 



SUBROUTINE FASTERCOUMT 

INPUT A polynomial pGR[3:;i] with exactly m>3 monomial terms, and a, fesR with a<b. 
OUTPUT The number of roots of p in the open interval (o, b). 
DESCRIPTION 

Step If m > 3 then use Sturm-Habicht sequences (as in, say, ILMOlp to count the number of roots in (a, b) and 
STOP. 

Step 1 Otherwise, let po '-—P, Pi '■=p' , let —p2 be the remainder of po/pi, and set i:=2. 
Step 2 If piGR then set K — i and GOTO Step 5. 

Step 3 Write Pi{x) —uix"^ — woi"", Pi-i{x) ■.= vix'''^ — vox''" (where ai >ao and 61 >6o), and let Cj :— ["^^^^j for 
iG{0,l}. 

Step 4 Replace i by i + 1, define pi{x) := -vi{^y^ x''^-"^^"^-"'''^ + ^o(Hi)=02.6o-<:o{ai-<io) ^ a^d GOTO Step 2. 
Step 5 Using recursive squaring (cf. RemarkUJ, evaluate and RETURN Na — Nb (cf. Definition|31l where A (resp. 
B) is {po{a), . . ■,PK{a)) (resp. (po(b), . . . ,px(&))). 

4. Correctness and Complexity: Proving Theorems [H El and El 
Proof of Theorem [H 

For open intervals, the second assertion follows immediately from Lemma ^ and Theorem El 
The case of closed intervals follows almost identically, save for an additional evaluation of 
/ at the end-points of the interval which increases the complexity bound by (a negligible) 
O(logD). 

The first assertion follows immediately from Theorems El and El and the second assertion 
of Theorem n which we've just proved. ■ 

Proof of Theorem El 

Henceforth, for convenience, we will say "time" in place of "arithmetic complexity." Clearly, 
it suffices to show that MNOMIALSOLVE is correct and satisfies the stated complexity bound 
when we set a* = a{D, m). Toward this end, let us first observe that Li{f) has exactly m — 1 
monomial terms and ^2(/) has <m — 1 monomial terms. Descartes' Rule then implies that 
Li{f) and L2{f) each have no more than m — 2 positive roots. Note also that /, /', and 
/" can all be evaluated using just 0(m log D) arithmetic operations, thanks to Remark [71 
These observations will be used implicitly throughout our proof. 

The correctness of MNOMIALSOLVE is then straightforward for m < 2. Also, for m = 2 
and -D>3, the complexity of MNOMIALSOLVE is clearly the same (asymptotically) as that of 
subroutine HYBRID with a*:=^^^, thanks to Proposition HI So we can assume m>3. 

In the special case where D = 2 and m = 3, note that Step 6 is but an implementation of 
completing the square. So, provided one shifts one's ^-approximations by 7p- while running 
the remainder of MNOMIALSOLVE, we can clearly attain the stated complexity bound for D = 2. 
So we can assume additionally that D>3 and now focus on Steps 10-18. 

In the special case 02 = 1, note that the change of variables x ^— > 1/x maps the interval 
(r — 6,r + 6) to ^21^2 (^ — S,r + 6) for r>6>0. An elementary calculation then reveals that 
5 <max {1, e'^/2} =^ for any 5-approximation y of x>e, l/y is an ^-approximation of 1/x. 
Note also that x"^/(l/x) = 03 + C2x"'''~°-^ + Cix"-'-* (since ai = after Step 4), and 03 — 02 > 1 
since 02 = 1 and a^ = D>3. So, provided one takes a reciprocal after running MNOMIALSOLVE, 
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we see that the worst-case time complexity of any instance with 02 = 1 is asymptotically no 
worse than that of the instances with 02 > 1. So we can assume additionally that 02 > 1- 

Now, by assumption, S{f) is dampened, as are all the inputs going into recursive calls 
of MNOMIALSOLVE. Note also that by our dampening assumption (and Rolle's Theorem), 
the intervals from Steps 17(b) and 17(d) indeed contain exactly one root of /". So our 
algorithm is well-defined and indeed finds ^-approximations to all the roots of / in [0,i?], 
assuming subroutines FASTERCOUNT and HYBRID are correct. The correctness of the latter 
two subroutines is proved below, so let us now concentrate on proving our main complexity 
bound. 

So let C{m) denote the time needed for MNOMIALSOLVE to execute completely for an input 
consisting of an / G ^(-D, m) (satisfying the dampening assumptions of our current theorem), 
an interval length of < R, and a precision e. Clearly, assuming subroutine HYBRID runs in 
time 0(mlog{D) log (a* log j)) (which is covered in a separate proof below), we then have 
that C{m) must satisfy the following recurrence relation: 



where the first term corresponds to finding the e-approximations of the roots of /' and /" 
(Steps 10, 13, and 14), the last term corresponds to the application of algorithm FASTERCOUNT 
to the intervals about these roots (Steps 13-16), and the O term corresponds to the appli- 
cation of algorithm HYBRID to the intervals between these roots (Steps 17-18). (Note that 
we have implicitly used the fact that /C(Z),m) and a{D,m) are non-decreasing functions of 
D and m.) In particular, for m>4, the last term of our recurrence is justified by Remark IHl 
Noting that we can regard R and e as constants, it is then clear that to find the asymptotics 
of our recurrence, it suffices to find the asymptotics of the recurrence 



A simple Maple calculation (rsolve(c(ni)=2*c(in-l)+a*m*(m-l)+b(ni) * (2*m-4) , c (m) ) ;) 
then shows that any such must satisfy Cm = 0(2"*(c2 + A) + ^2^=1^^"^ J ^U))- So, 
knowing that C(2) = 0(log(D) loglog j) |Ye94| Pg. 279], we must then have that 



Proof of Theorem ^ 

The complexity bound follows immediately from the upper bound on K, assuming algorithm 
FASTERCOUNT is correct. The bound on K and the correctness of algorithm FASTERCOUNT are 
proved below. ■ 

At this point, we need only show that the subroutines HYBRID and FASTERCOUNT are 
correct, and that FASTERCOUNT satisfies a sufficiently good complexity bound. 

Proof of Correctness and Complexity Analysis of HYBRID: 

First note that by assumption, there is a unique j G {0, . . . ,2^~^} such that {jS^e, P^^^e] 
contains the sole root of in (0,i?). It then becomes clear that Steps 0-2 are merely an 
implementation of bisection that finds this j: In particular. Steps 0-2 essentially find the 
binary expansion of j (from most significant bit to least significant bit), and x = (3^ once 
Step 3 is reached. By Theorem HI x is then an approximate root of 0, and Step 3 indeed 
finds an e- approximation of our root. 



C{m) < 2C{m — 1) + O I m{m log D) log I a* log 




Cm<2cm_i + TV? A + (2m - A)B{m). 



C{m) = 0[2"'{ \og{D) log a* log - + ni^K,{D, ni) 
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As for the complexity of our algorithm, it is clear that M= \\0g2 (log2 7) ~ log2 log2 /3] and 
we thus need only O (log (log j) — log log /?) evaluations of (p (and inequality checks) until 
Step 3. Step 3 then clearly takes just O (log log j) arithmetic operations and evaluations of 
(f) and 0'. 

To conclude, recall the elementary inequalities ^^^^^^^ < x, valid for all a; > 0. Observe 
then that — log log (3 = log < log a*. Since an m-nomial of degree D implies that and 
0' can be evaluated within 0(m log D) arithmetic operations (cf. Remark [71), we are done. 
■ 

Proof of Correctness and Complexity Analysis of FASTERCOUNT: 

The correctness of FASTERCOUNT when m > 3 is clear from any of the standard references 
|Roy96[ ILMOlj so there is nothing to prove. 

For m = 3, correctness follows easily from Lemma Q upon noting a simple fact: The 
remainder of ^^^„^~^°^' ° is nothing more than the reduction of fix^^ — fnX^° modulo Utx"'^ — 
MoX°°. In particular. Step 2 of FASTERCOUNT is just a compact representation of a sufficient 
number of applications of the identity x"'^ = ^x"" to reduce —{vix^^ —VqX^^) modulo Uix""^ — 
uox°'° to a polynomial of degree <ai. So correctness follows immediately. 

As for the complexity bound, we need only observe that (by recursive squaring) every 
execution of Step 2 and Step 3 takes only 0{logD) arithmetic operations, and that K = 
0{logD) at the termination of the algorithm, i.e., there are only 0{logD) remainders in our 
Sturm sequence. The first assertion is clear, so to prove the latter assertion we will need to 
prove a technical bound on the exponents which occur in our remainder sequence. 

In particular, let ii be the absolute value of the difference exponents of pi (we set : = if 
Pi is monomial). Note that pi a monomial =^ pi^i is a monomial, and thus pi+2 is constant. 
Note also that, via the definition of long division, we have £j+2 < — ^il? with equality 
occurring iff pj_(_2 is a binomial. So if we can show 

(*) mm{ii+2,ii+3}<£i/'2 

for all i, then we will easily obtain our bound on K (since io,ii < D and all the £j are 
integers) . 

To prove {-k) observe that < =^ > (if 7^ 0) and thus £^+3 < 

Similarly, ii+i>jii =^ ^i+2<|^j- So < 3 [logg -D] + 2 = 0(logD) and we are done. ■ 

5. The Proof of Theorem [21 

Dividing by a suitable monomial term, we can clearly assume without loss of generality 
that / has a constant term . Note then that C is diffeomorphic to a line iff / is not the 
square of binomial, via |LRWn3[ Prop. 2 and Lemma 1]. So there are actually no isolated 
singularities for trinomial curves. 

It is then clear that X IS cL vertical tangent of C =^ X IS cL root of := (/, ^). The roots 
of Fy with X or y coordinate equal to or i? can be found easily by suitably fixing one of 
the variables and applying Theorem^ so let us concentrate on the roots of in (0, R)"^. 

By |RojOOa[ Lemma 4.2] and vdKOOj we can then find (for any constant 5 > 0) a mono- 
mial change of variables, using just 0(log^''''' D) bit operations and O (\og{D) log [D log j)) 
arithmetic operations, such that the resulting sparse polynomial system G := {gi,g2) has 
Qi gM[xi, X2] a trinomial, (72 a binomial, and all underlying exponents have bit-length 

0{\og^D). We can then solve g2 to accuracy min < 1 [ via Theorem [TJ back-substitute 
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the roots into gi, solve the univariate speciahzed gi to accuracy min je'^^'^^'', l|, then invert 
our monomial map via IRojOOat Main Theorem 1.3] and |vdKnnj to recover the roots of F„, 



using just 0{\og^~^^ D) additional bit operations. So we can indeed find isolated vertical 
tangents within the stated complexity bound. 

To conclude, note that |LRW03| Lemma 9] tells us that if x is an isolated infiection point 
of C then [dff- {82 ff - 28182/ ■8if-82f + 8^f-{8iff]=0, where 8i:=£-. In the special case 
of a trinomial, say f{xi,X2) := 1 + Ax1x2 + -Bx^Xj, the preceding polynomial in derivatives 
is exactly: 

-ab (a + b)R^ + {l?c^ + a^(f -b^c-2 abd -a^d-2 abc - 2 abed) R^S 

+ {b^c^ - 2 abed -ad^ -2 acd -2bcd- bc^ + a^d^) RS^ -cd{c + d) S^, 

where R:=Ax'^X2 and S := Bx\x2-, thanks to a Maple calculation. (Note that we must also 
have 1 + i? + S' = since we're assuming / = 0.) So we efficiently solve for x as follows: (1) 
solve the above cubic for ^ up to accuracy min \^e^^^\ l} via the main algorithm of jNR96j 



or [BP94] . (2) solve for (i?, S) via the additional relation l + i? + S' = 0, (3) solve the resulting 
binomial systems in (i?, S,xi,X2) for x via [RojOOat Main Theorem 1.3] and |vdK00j . So we 
can e-approximate the infiection points of C within the stated asymptotic complexity bound 
as well, and we are done. ■ 



6. The Proof of Theorem [HI 

Definition 8. For any aeR and /cGNU {0}, let {a)k: = a{a-l) ■■■{a-k + 1). (So (a)o = 1, 
{a)i = a, oGN => (a)a = a.!, and a, A;GN with k>a =^ {a)k = 0.) o 

Proof of Theorem HI 

The first bound is immediate since J^(-D, m) C JF(Z}, m + 1). 

The formula for a{D,2) follows easily since m = 2 and /(O) 7^ implies that /' has no 
positive roots. This in turn implies that N{f) = 1, so the only interval we (may) need 
to majorize over is (mo(/), Mi(/)) = (0,cxd). In particular, the quantity we majorize is just 
X7(/, x), and the formula for a{D, 2) then follows immediately from Proposition^ So we can 
assume m = 3. Since 7(/, x) = 7(0/, x) for any nonzero constant c, we can also clearly assume 
that f{x) = x""^ — Ax""^ + B where A and B are nonzero real constants and 03 = > 02 > 1. 
Since we assume /'(0)=0 in the definition of a{D,m), we must also have D>3 and 02 >1. 

Next, note that A < =^ /' has no positive roots. So, similar to the m = 2 case, one 
either majorizes over the empty set or (0, 00). In the former case there is no contribution to 
a{D, m), while in the second case, we see that we are again majorizing X7(/, x) over (0, 00). 
Since the latter supremum is no greater than by Propositions El and |3J we can can now 
assume additionally that A>0. 

Now note that 

fix) = asx^'-^ - a2Ax^^-^ 

= asX-'-'(^x'''--'~^y 
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and 



Clearly then, Xi :- 



f{x) = 03(03 02^(02 

03(03 «2A(o2-i: 

\ 1/(03-02) 

a2A \ 



03(03 - 1) 

is the unique positive root of /' and X2 



as J ...^ ^...^^^ t-^^...,^ . .J y asias-l) J 

is the unique positive root of /". Furthermore, since 1 < 02 < 03, it is clear that < 0:2 < a;i. In 
particular, we see that we are left with a collection of 3 intervals to consider in the definition 
of 0({D, m): 

{(0,a;2), (2:2, 2^1), (a;i,oo)} 
and the quantities we majorize on these intervals are, respectively: 

{xi{f, x), {xi - x)-f{f, x),{x- Xi)-f{f, x)} . 

The intervals (x2,xi) and {xi,oo) can be handled via a unified calculation, so let us 
complete our analysis by examining the intervals {x2,oo) and (0, a;2). 

The Interval (xa, 00): 

Note that f{y + Xi) = {y + Xi)"'' - A{y + xi)"^ + B and thus 

/W(y + xi) = {aMy + ^ir'-'-{(^2)kA{y + x,r-\ 
So, letting 0(|/) : = f{y + Xi), we have for all k>2: 



(1) 
(2) 

(3) 



iy) 



i/(fe-i) 



k\{as{y + xiYi-^ - a2A{y + a;i)"2-i) 



(03)/= 



i/(fc-i) 



i/(fc-i) 



fc! 



03 



i/(fc-i) 



{y + xiY-^ ((y + a:i)"3-"2 



i/(fc-i) 



which in turn is bounded above by 



(4) 



since 



\y\ ' ' 



1 



(<>2)kA 

("Yk 



{y + xi 



\k-l 



{y + xi)"3-"2 - xl 



,03—02 



i/(fc-i) 



(03)fc 
fe!o3 



l/(fc-l) 



< 03-1 
— 2 • 



Now let z:—y/xi and let p be ■ ^^^^^ 
for y^—xi (or z^—\, equivalently) we have: 



^22dL- \ / xi or according as /c < 02 or not. Then 



(5) 

(6) 
(7) 



(i/ + xi)'=-i (j/ + a;i)<'3-«2_xJ3-<'2 



i/(fc-i) 



(^ + l)'=-i 



+ l)«3-a2 _ 1 



l/(fc-l) 



{Z + l)«3-a2 _ pa3-a2 



l/(fc-l) 



(«+l)'=-l l + («+l)H h («+ l)"3-«2-l 

{fc-2)/(fe-l) _|_ ^^(,3-02 _ pa3-02 



Z+1 



:+l) + (2+l)2--- + (z+l)" 



l/(fe-l) 



where the last equality follows from the identity = = i+{z+i)+-\{z+i)'^ --^ ■ 
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Since = 1 — is clearly an increasing function of z for all z > — 1 (or y > —Xi, 
equivalently), and since 0>— ^>— 1, we must have for aA\ y > X2 — Xi (or 2; > ^ — 1, 
equivalently) that 



z+l 



(fc_2)/(fc-i) r (fc_i)/(fc_2) 

< max < I 1 1 , 1 

X2 



, 11 > N (fe-l)/(fe-2) 



a2 — 1 



- 1 



(*) 



< max < 1,1 

02 — 1 



Furthermore, since "^^^ < 1 for 2 < < ao — 2, we must have — > p and thus the last factor 
of Equality ((7j) is bounded above by 1. So by Equalities (H}-© and (^) we obtain that 



0W(y) 



i/{k-i) 



< - max{(a3 - l){as - 2), - 1} < {D - 1){D - 2)/2, 



since D >3. Thus, the contribution of the interval {x2,oo) to a{D,m) is no more than 
(D-l)(D-2)/2. ■ 



The Interval (0, ^^2): Similar to the last case, we have for all a; G (0,^2) that 



X) 



klf'ix) 



l/(fc-l) 



ia3)kX^^ - {a2)kAx' 



a2—k 



(a3)fca;'^3-"2 - (as)^^ 



i/(fc-i) 



i/(k-i) 



1 


(as)^ 


i/(fc-i) 


^03—02 


X 


A;!a3 







l/{k-l] 



< 



1 03-1 

X ' 2 



'1 

ag-aa _ (a2)fcA VC^-l) 



Xa3-a2 _ a;^ 



(13—02 
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So, for 0<z<X2/xi = (f|5i)^/("3-a2)<i^ h^ve 



X 



03—02 



(a3)fc 



^03—02 



X 



03—02 



l/(fc-l) 



^03—02 



^03— 02 



^03—02 



l/(fe-l) 



l/(fc-l) 



< max 



< max ■ 



02 





1 


_ p03-a2 


l/(fc- 




1 


^03—02 




"j 






1 


l/(fc- 


"1 


1 


^03— <I2 







(recalling that p<l) 
(recalling that p< 1) 



< 

< a2 - 1 

< 03-2. 



i/(fe-i) 



Thus 



/W(x) 



A;!/'(x) 



i/(fc-i) 



< 



1 (a3-l)(«2-l) 



y 2 

and the quantity we need to majorize on this interval is no more than {D — !)(£) — 2)/2. ■ 
Putting together the bounds we've found over our two preceding intervals, we see that our 
upper bound for a{D, 3) holds and we are done. ■ 
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